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1. Introduction 



Hadron colliders, such as the Fermilab Tevatron and the CERN Large Hadron Collider (LHC), 
bring us to the high-energy frontier for new discoveries of physics beyond the Standard Model 
(SM). Yet they are QCD machines: The strong interaction of the colored partons, quarks and 
gluons, leads to overwhelmingly dominant event rates in hadronic collisions. Thus any new 
particles participating in the QCD interaction will be strongly produced with favorable rates. 

It is often said that the LHC is a "gluon factory," as initial states involving gluons are 
the major contributors to the majority of the events, particularly for relatively light final 
states. On the other hand, due to the abundance of valence quarks in the proton at large 
momentum fraction x, production via initial-state quarks becomes more important for rather 
heavy particles. Since the LHC is a pp collider, production of heavy particles in qq collisions 
will be significantly enhanced by the valence-valence component of the initial state. A simple 
group theory consideration indicates that quark-quark scattering can produce bosonic color 
states 6© 3. Particles in these color representations would be rather exotic states beyond the 
SM spectrum. 

Incidentally, there are indeed theoretical motivations to consider those exotic color states. 
Examples include scalar quarks in supersymmetric (SUSY) theories with iZ-parity violation 

and color-sextet scalar diquarks Q.^ Diquark states of exotic quantum numbers have 
recently been classified in Ref. [^. Resonant production of color-antitriplet scalars and vectors 
at the LHC have been studied in Refs. [Q, ^, ^. Pair production of color-sextet scalars at 
the LHC has been studied in Ref. [0]. In this paper we consider the single production of 
these colored scalar states at hadron colliders including the next-to-leading order (NLO) 
QCD corrections. We also present their transverse momentum spectrum by including the 
leading-log soft gluon resummation to all orders in the strong coupling. These calculations 
have been done previously for single scalar top production in i2-parity violating SUSY |^] ; we 
have expanded upon that work. 

2. Model And Constraints 

In this paper we are not concerned about the electroweak structure of the diquark, but for 
completeness Table || lists the possible electroweak quantum numbers of the diquark and 
its coupling to standard model quarks. Q is the SU{2)l quark doublet, and U(D) is the 
up(down)-type SU{2)l quark singlet. The diquark may be either an antitriplet or a sextet 
under SU{3)c, independently of the electroweak quantum numbers. Also, the diquark may 
be either a scalar or vector boson. For simplicity and due to theoretical motivations j^, ^, 
we concentrate on the scalar case. 

Any states with nontrivial SU{3)c color quantum numbers will interact with gluons 
and can thus be produced via gauge interactions. However, to understand the properties of 

^In this work, the word "diquark" refers not to a bound state of two quarks, but rather to a fundamental 
boson carrying the quantum numbers of two quarks. 
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Table 1: Electroweak quantum numbers of possible diquark states and their allowed couplings to SM 
quarks. In each case the SU{3)c quantum number may be either 3 or 6. Under SU{3)c x SU{2)l x 
U{1)y, Q has the quantum numbers (3,2,1/6), U has the quantum numbers (3,1,2/3), and D has 
the quantum numbers (3, 1, —1/3). 
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exotic colored states, we must study their model-dependent aspects. We consider a model 
parameterization in which, after electroweak symmetry breaking, a colored scalar diquark 
(denoted by D henceforth) couples to quarks via the Lagrangian 



c = 2V2 



KrD'qa[\LPL + \RPR)q\; , (2.1 



where Pl,r = (1 T 75)/2 are the left- and right-chirality projection operators, is the 
conjugate quark field, and the sum over quark flavors has been suppressed. The scalar 
transforms according to either the sextet or antitriplet representation of S'f7(3)c, and the 
Ki"'^ are Clebsch-Gordan coefficients coupling this representation to two triplets (see Ap- 
pendix^ for more details.) Also, for different flavor combinations the couplings \l,r may be 
completely independent, or they may be constrained by some underlying model. Note that 
in the antitriplet case the couplings must be antisymmetric in flavor. 

Existing data impose constraints on the couplings A in specific models. For a review of 
constraints on i2-parity-violating SUSY models see Ref. In the case of a sextet diquark 
with a mass of a few hundred GeV to 1 TeV, coupling to right-handed up-type quarks, 
constraints arising from D^-D^ mixing and the non-strange D meson decay to tttt require |^ 

A^", A^^<0.1andA^^«0. (2.2) 

Additionally, since the left-handed Cabibbo-Kobayashi-Maskawa (CKM) matrix is known, 
the left-handed couplings A^, are tightly constrained due to minimal flavor violation (MFV). 
However, these constraints do not apply to the right-handed couplings \r, and our results 
are sensitive only to the combination A^ = A| + A|j. Since we do not assume any particular 
model, we do not concern ourselves with these constraints except to observe that Yukawa 
couplings of order 0.1 are allowed for the diquark masses of our interest. 

3. Leading-Order Scalar Production 

At leading order a scalar diquark is produced via the process 

q{pi) + q{p2)^ Dil), (3.1) 
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where pi, P2, and I are the momenta of the particles. The invariant amphtude for this process 
is 

= -2V2icl)iU^''{pi)K\b{XLPL + XrPb)C^u\p2), (3.2) 

where C is the charge conjugation matrix and (pi is the color wavefunction of the diquark. 
The Born-level partonic cross section is then 

^BoUs) = ^^^5(1 - r) = ^6(1 - r), (3.3) 

where s = {pi +^2)^, = mj^/s, mo is the mass of the diquark, Nq = 3 is the number 
of colors in the fundamental representation, and No is the dimension of the diquark's color 
representation. That is, Nd = 6 for a sextet diquark and = 3 for an antitriplet diquark. 

The hadronic cross section is given by a convolution of the partonic cross section with 
parton distribution functions (pdfs): 

f^LO = ^2 (g (3.4) 

c 

where S is the hadron collider's center-of-mass energy squared, tq = mj^/S, q{x) is the quark 
distribution function (with a flavor sum again suppressed), and (8) denotes the convolution, 
defined by 

{h®f2){x)= / dxi / dX2 5{xiX2- x)fi{xi)f2{x2), (3.5) 



which is the /1/2 parton luminosity. 

4. Next-to-leading Order QCD Corrections 

Next-to-leading order QCD corrections to scalar diquark production arise through virtual 
gluon loops, the real-gluon emission process qq — )• gD^ and the gluon- initiated process gq — )• 
qD. We give a general outline of the calculation in this section, and provide more details in 
Appendix 

We do not consider the decay of the diquark. Since the diquark carries color charge, 
a calculation taking its decay into account would include QCD corrections involving the 
exchange of color between the initial and final states. Thus the NLO QCD corrections do 
not factorize perfectly into production and decay. However, the effects of this imperfect 
factorization are suppressed by F^i/m^), where T^) is the decay width of the diquark. In the 
present work we neglect such effects, and treat the diquark as a stable particle produced on 
shell. 

4.1 Virtual corrections 

The one-loop diagrams contributing to the process qq ^ D are shown in Fig. These virtual- 
gluon loops give rise to ultraviolet (UV), soft, and collinear divergences. All the divergences 
are regulated using dimensional regularization in 4 — 2e dimensions. 
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Figure 1: Feynman diagrams for virtual gluon corrections to qq ^ D. 



The UV divergences of the loop diagrams necessitate renormalization of the quark and 
diquark wavefunctions and of the qqD vertex. We perform this renormalization in the MS 
scheme. The wavefunction renormalization constants are found to be 

^2 - ^ " An r(l-e)7' ^2 - 1 + 27r r(l-e)7' ^^'^^ 
and the vertex renormalization constant is 

In these equations, Cp and Cd denote the eigenvalues of the quadratic Casimir operator of 
SU{3) acting on the fundamental representation and on the diquark representation, respec- 
tively. For a sextet diquark, Cd = 10/3, and for an antitriplet diquark, Cd = 4/3 = Cp. 
Since we are considering the single production of an on-shell diquark, the renormalization of 
the diquark mass is not needed for the calculation. 

The effect of the self-energy diagrams and counterterms is to replace the coupling in the 
Born cross section by 

A^ZA(Z|)-^(Zf)-i/2A, (4.3) 

where A on the right-hand side denotes the renormalized coupling. Once this replacement is 
performed, the remaining virtual corrections are given by the triangle diagrams in Fig. |l|. 

For the complete NLO QCD treatment, we also include the running coupling A(^|j), 
given by 

where fiR is the renormalization scale. 
4.2 Real gluon emission 

We must also consider the NLO corrections due to the radiation of an additional parton 
into the final state. The process qq — >• gD, shown in Fig. features both soft and collinear 
divergences. The soft divergences cancel against those from the virtual correction, and the 
collinear divergences are absorbed into the pdfs using the MS factorization scheme. 
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Figure 2: Feynman diagrams for qq gD 



The total hadronic cross section corresponding to the partonic subprocess — t- D + X 
is then given at NLO by 



a 



m 

NLO 



27rA2 Nd dr 



a. 



T \ T 



5(1 -r) 



m 



+ 2Cf 
-Cd 



2(l + r2)fi^^^i^l + 



1 + 
(1-r). 



+ 



1 - r 
27r2 



vr 



3 , m 



3 2ln9^-l)5(l-T) + l-r 



(4.5) 



1 5(1 - r) 



where /^i;' is the factorization scale, and Pqq{T) = Cf\{1 + t)^/(1 — t)]^ is the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting function B. 



4.3 Gluon-initiated process 




Figure 3: Feynman diagrams for gq — > qD. 

The 55 partonic channel contributes for the first time at NLO, via the process gq — )■ qD^ 
shown in Fig. ^. This process has only collinear divergences, which again are treated via 
the MS factorization scheme. The hadronic cross section corresponding to the gq partonic 
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channel is 



91 

NLO 



[ —{9®q + q®g){— 



X^P,g(T)ln ^ ^ 



i(l-r)(3-7.) + §2 



TlnT + -(l-T)(l+2r) 



(4.6) 



where Pqgij) = | ["^^ + (1 - t)^] is the DGLAP sphtting function. 

The total cross section for inclusive scalar diquark production is given by the sum of 
Eqs. Q and 





Figure 4: Results for -pp collisions at a center of mass energy of 2 TeV (a,b) and -p-p collisions at 
10 TeV (c,d) with various initial states. The total leading-order (dot-dash) and next-to-leading order 
(solid) cross sections are shown for both (a,c) the antitriplet and (b,d) the sextet diquarks. For all the 
above plots the factorization scale, renormalization scale, and diquark mass are set equal. 

5. Numerical Results 

In this section we present numerical results for the total cross section for production of a 
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scalar diquark, for various diquark masses, center-of-mass energies, and initial states. We 
use the CTEQ6L set of parton distribution functions |10] for leading-order results, and the 
CTEQ6.1M set HI for NLO results. Motivated by i2-parity violating SUSY |], for the 
antitriplet diquark we only consider initial states that can result in a final state with charge 
—2/3 or +1/3. Also, following a partially unified Pati-Salam model for the sextet case 
we only consider initial states that can result in a charge +4/3 final state. Note that in the 
antitriplet case there are no identical initial partons since the couplings must be antisymmetric 
in flavor. Unless otherwise noted, the factorization scale, renormalization scale, and diquark 
mass are set equal. 

Due to the much larger parton luminosity for the valence quark initial states at the LHC, 
we only consider qq contributions. The inclusion of the qq initial state would double the 
diquark production cross section for all initial states at the Tevatron (a pp collider) and for 
pure sea-quark initial states at the LHC. If a valence quark is in the initial state at the LHC, 
the qq initial state contribution is much smaller. For example, at the LHC at 14 TeV with 
a diquark mass between 500 GeV and 1.5 TeV, the uu initial state has a cross section 10-80 
times larger than the uu initial state. 

The leading-order and NLO total cross section results for the antitriplet (sextet) diquark 
are shown in Fig. |4(a)| (Fig. [4(b)D for the Tevatron at 2 TeV, Fig. |4(c)| (Fig. f4(d)D for the 
LHC at 10 TeV, and Fig. (Fig. ^(b)l) for the LHC at 14 TeV. We have factored out 
the coupling for a model-independent presentation. At lower diquark masses and not too 
small coupling, the Tevatron would produce scalar diquarks at reasonable rates, reaching a 
cross section of the order of picobarns near m/j ~ 700 GeV. At the LHC at 10 TeV the 
production cross section drastically increases over the Tevatron by a factor of 10^-10^ for 

= 500-1000 GeV. The production cross section at the LHC at 14 TeV increases slightly 
over the 10 TeV cross section. Also, the NLO cross section increases over the leading-order 
cross section in all cases. If the couplings between the quarks and diquarks are not too small, 
the diquark will be produced at favorable rates. 

One generally expects the cross section for production of a sextet diquark to be larger 
than that of an antitriplet, due to a larger color factor in Eq. (3.4). Thus, for example, the 
cross section of the sextet production from uc is larger than that of the antitriplet production 
from ds. The comparison of uu sextet production to du antitriplet production is an interesting 
case. Apart from the larger color factor, uu receives a relative enhancement due to the larger 
up-quark pdf. However, du also receives an enhancement, due to the combinatorics of the 
initial state: the du luminosity is d®u + u®d, while the uu luminosity is simply u®u. These 
luminosity enhancements approximately offset one another, so that the uu cross section is 
approximately twice that of du. 

An important quantity for NLO calculations is the i('-factor, defined as the ratio of the 
NLO cross section to the leading-order cross section. Figures 5(c), 5(d) show the X-factor 
as a function of the diquark mass for the antitriplet and sextet cases, respectively. In the 
antitriplet case the ii'-factor is between 1.31 and 1.35. For the sextet case, the X-factor 
runs from around ~ 1.27 to ~ 1.32 for the up-charm initial state, is around 1.22 for the 
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Figure 5: Results for pp collisions at a center of mass energy of 14 TeV with various initial states. 
The total leading-order (dot-dash) and next-to-leading order (solid) cross sections are shown for both 
(a) the antitriplet and (b) the sextet diquarks with various initial states. Also shown are the if -factors 
for (c) the antitriplet and (d) the sextet diquarks with various initial states. For all the above plots 
the factorization scale, renormalization scale, and diquark mass are set equal. 



up-up initial state, and is between ~ 1.49 and ~ 2.15 for the charm-charm initial state. The 
purely sea-quark initial states (cc) show more dependence of the iC-factor on the diquark mass 
because of the differences in the factorization-scale dependence between the leading-order and 
NLO pdfs. For initial states involving up and down quarks, the iT-factor for the sextet case 
is generally smaller than that of the antitriplet case. This surprising result is due to a partial 
cancellation between the two color structures [the Cp and Cd terms in Eq. ( [1.5[ )]. 

The factorization and renormalization scale dependence is shown in Figs. |6(a) , 6(b) for a 
500 GeV and a 1 TeV sextet diquark, respectively, and in Figs. 3(c) , 6(d) for a 500 GeV and a 1 
TeV antitriplet diquark, respectively. For all the figures, the factorization and renormalization 
scales were set equal and varied from m£)/4 to m/). As expected, the scale dependence of the 
NLO cross section is less than that of the leading-order cross section for all cases. 

We have also computed the rapidity distribution of a 1 TeV sextet diquark at the LHC. 
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Figure 6: Leading-order (red) and next-to-leading order (black) factorization and renormalization 
scale dependence for pp collisions at a center of mass energy of 14 TeV and various initial states. 
Figures (a) and (b) show the factorization and renormalization scale dependence for a 500 GeV and 
1 TeV sextet diquark, respectively, and figures (c) and (d) show the scale depence of a 500 GeV and 
1 TeV antitriplet diquark, respectively. The factorization scale and renormalization scales are set equal 
in all plots. 



The shape of the NLO distribution is virtually identical to that of the leading-order distribu- 
tion; the only significant difference is in the overall normalization. Thus, for purposes of the 
rapidity distribution, the effect of the NLO corrections is well described by a simple /C-factor. 

In Table | we tabulate the cross sections at leading order, leading order with NLO pdfs, 
and full NLO, as well as the exclusive cross sections for production of a diquark with an extra 
parton (gluon or antiquark) in the final state. This final-state parton is required to have a 
transverse momentum of at least 20 GeV. When a parton appears in the final state, it is more 
often a gluon than an antiquark, by a ratio of about 3 : 1 if = 500 GeV, or about 8 : 1 if 
m£, = \ TeV. This is due to the increasing dominance of the valence-quark parton luminosity 
over the gluon parton luminosity as the diquark mass increases. 

As for the possible observation of the diquark signal, the decay of the diquark to two 
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Table 2: 14 TeV LHC total cross section (in nanobarns) at leading order, leading order with NLO pdfs, 
full NLO, diquark plus a gluon, and diquark plus an antiquark. A minimum transverse momentum of 
20 GeV is required for final-state partons. For the sextet (antitriplet) diquark case the initial state is 
uu {du). In this table the coupling A is taken to be 1; the cross sections scale as A^. 
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light-quark jets may be masked by SM QCD dijet events. The decay to two tops can be more 
distinctive, in particular for the like-sign charge combinations of electrons and muons from 
the tt decay, although the branching fractions of the diquark to two tops may not be too large 
and will be model-dependent 



6. Soft Gluon Resummation and the Transverse Momentum Distribution 



In Table |2|, we listed some cross section values for a diquark plus a final-state parton. These 
rates are very sensitive to the cutoff on the parton transverse momentum px- For transverse 
momenta above other scales in the process, say the diquark mass, fixed order perturbation 
theory in as gives reliable results. If the transverse momentum is much smaller than the 
diquark mass, the series in ag is replaced by a series in ag 

— ^ the expansion becomes non-perturbative, but it is possible to resum the terms that 
are at least as divergent as ■ Using the results in the previous section, we can obtain the 
full pt distribution of the diquark at the leading log to all orders in as ■ 

For the production of strongly interacting particles, the formalism for the resummation 
of small transverse momentum logarithms has not been fully developed |jl^. Although it has 
not been proven, we assume that the final-state and initial-state radiation can be resummed 
in the same way. This assumption has been made previously for tt production at the Tevatron 
[ p^ ] and scalar top quark production at hadron colliders Q. 

Following the usual procedure fl^ , the leading-order asymptotic cross section in the low 
Pt limit is 



a. 



dpT dy ^ -K SpT 



2 

Pt 



(6.1) 
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where y is the diquark rapidity, 2 = {mD/\^)e^^ are the dominant terms at low px, and 
the coefficients A, B are found to be A = 2Cf and B = —{'iCp + Cd)- Our results for A and 
B agree with the previous result for single top squark production [Q] . The resummed formula 
for diquark production is 



(Pa' 



dpxdy 
W{b) 



2PT 

S 



db-Jo{bpT)W{b), 



(6.2) 



exp 



2tt 



Aln- 



m 



D 



+ B 



where Jq{x) is a zeroth-order Bessel function, W{b) is a Sudakov form factor from the sum 
of large logs in impact parameter space, the canonical value of 60 is 2e~'^^ , and the pdfs are 
evaluated at a factorization scale oi = bo/b. If the running of ag is kept to leading order, 
the integral in Eq. ( |6.2[ ) can be evaluated. Also, Aqcd was set so that the leading-order 



as{M'^) attained the correct value. 

For b > 1/Aqcd confinement sets in and diverges. To cut off the divergence and 
parameterize the non-perturbative effects, we make the replacement IS, 20| 



W{b) W{b^)exp 



,2 1 "max"lD 

b g2 In 



6* 



(6.3) 



(2GeV) [21]. There are other methods 



where gi = 0.14 GeV^, g2 = 0.54 GeV^, and 6ma= 
of parameterizing the non-perturbative effects but the results did not significantly differ 
from this parameterization. 

The coefficient A = 2Cf is of the usual form found from resummation of the Drell-Yan 
production [19, 20 [, while the other coefficient B = —{'iCp + Cd) contains the usual term 
—^Cp found in the resummation of Drell-Yan production []l9| , |20| [ and a new term proportional 



to Cd- In order to investigate the sensitivity of the results to the assumption that initial- and 
final-state radiation resum in the same way, we have also computed the resummed distribution 
without final-state radiation, i.e. with B = —'iCp- We find that this causes the peak in the 
resummed distribution for the antitriplet du (sextet uu) initial state to decrease by about 
20% (45%). Also, the peak in the resummed distrbution for the 500 GeV (1 TeV) diquark 
shifts to the right by about 2 GeV (1 GeV). Since C^) is larger in the sextet case than in the 
antitriplet case, setting B = —SCp has a larger effect on the resummed distribution in the 
sextet case. 

The resummed distribution has only accounted for the terms in the perturbative ex- 
pansion that are at least as divergent as p^^ ■ In the high-pT region other terms become 
important and the resummed distribution is inaccurate. To match between the perturbative 
and non-perturbative regions the total px distribution can be defined as 



dpT dy dpT dy 



+ f{PT) 



dpT dy dpT dy 



Kpt) 



1 



1 + (pt/Pt'"*''^)^' 



(6.4) 
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where p^'^ is a scale above which the the perturbative distribution is accurate, and f{pT) 
is a matching function. The function / is needed, since at large pT the resummed and 
asymptotic distributions are invalid and their difference can be larger than the perturbative 
distribution's value. The large power of pT in the denominator of /{pt) is chosen so that the 
matching function quickly goes to zero as px increases above p™'^*'^'^, and quickly goes to one 
as pt goes to zero. 
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Figure 7: px distributions or diquarks with masses of 500 GeV and 1 TeV. (a) shows the pT distri- 
bution of an antitriplet diquark with du initial state, and (b) shows the pt distribution of a sextet 
diquark with uu initial state. Both (a) and (b) show the aymptotic (dashed), perturbative (dot-dash), 
resummed (dotted), and total (solid) distributions. 

Figure 1^ shows the pT distributions for both a 500 GeV and 1 TeV diquark. The antitriplet 
Pt distribution with down-up initial state is shown in Fig. 7(a) , and the sextet pT distribution 



with up-up initial state is shown in Fig. 7(b) . For all the distributions the typical value of 



pmatch — jyj^^3 T^g^g ugg(j, At low pT the pcrturbative and asymptotic distributions cancel 
and the total distribution matches the resummed distribution, while at high pT the matching 
function goes to zero and, beginning around px = iud/S, the total distribution converges to 
the perturbative distribution. As is typical, the low px distributions peak around 5-8 GeV. 
Hence, the pT distribution no longer diverges as px ^ and the perturbative distribution 
becomes reliable around m£,/2>. 



7. Summary 

The LHC is a hadron machine, hence any new particle participating in QCD interactions can 
be produced with favorable rates. Due to the gluon's high parton luminosity, the LHC is typi- 
cally refered to as a "gluon factory." For rather heavy final states the gluon parton luminosity 
decreases relative to the valence-quark parton luminosity, hence valence-quark scattering is 
still important. Quark-quark annihilation can result in colored sextet and antitriplet scalars, 
so-called diquarks. 
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We calculated the NLO corrections to the single production of these exotic colored states. 
The production cross section at the LHC was favorable, provided that the coupling between 
the diquark and SM quarks is not too suppressed. For most initial states and a diquark mass 
between 500 GeV and 1.5 TeV, the K-factor for the antitriplet diquark is between ~ 1.31 
and ~ 1.35, and for the sextet diquark the if-factor is around 1.22 for the up-up initial 
state and runs from ~ 1.27 to ~ 1.32 for the up-charm initial state. If the initial state is 
composed exclusively of sea quarks, the i^- factors are much larger due to the differences 
between the leading-order and NLO pdfs. Also, the NLO calculation was found to decrease 
the renormalization and factorization scale dependence for all initial states and both the 
antitriplet and sextet cases. 

The soft gluon resummation for low transverse momentum was also calculated for the 
antitriplet diquark from the down-up initial state and the sextet diquark from the up-up 
initial state with diquark masses of 500 GeV and 1 TeV. At low px the total pT distribution 
matches the resummed distribution, peaking around 5-8 GeV. For pT ^ ra^i/'i the total 
distribtution converges to the perturbative distribution, indicating that for pT above m£)/3 
the perturbative distribution is reliable. 
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A. Color 

The product of two fundamental representations of SU{3) is 3 <S> 3 = 6 © 3, so a boson 
produced in qq fusion can be either a color sextet or an antitriplet. Since the sextet especially 
is rather unfamiliar, we include in this Appendix some technology for working out the color 
factors involved in our calculation. While the gauge group SU{3) is of primary interest, we 
present results for arbitrary SU{Nc)- 

A.l The fundamental representation 

We denote a vector transforming under the fundamental representation of SU{Nc) by a raised 
Latin index: m°. An object with a lowered index, Ua, belongs to the antifundamental. The 
generators of SU{Nc) in the fundamental representation are denoted by t'^^f,, or just in 
matrix notation. Uppercase Latin indices belong to the adjoint representation; they can be 
freely raised and lowered using the Kronecker symbols S^^ and Sab- The generators are 
given by = ^X^, where the A"^'s are the familiar Gell-Mann matrices, and they satisfy the 
commutation relation 

[i^,i^] = z/^^%, (A.1) 
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where Z"^^*^ are the structure constants of SU{Nc)- The quadratic Casimir operator is given 

by 



t^tA 



and the generators satisfy the orthogonahty relation 



2 



and the completeness relation 



1 



m 



(A.3) 



(A.4) 



A.2 Clebsch-Gordan coefHcents 

The generators T"^ in the product of two fundamental representations can be written as 

rj^Aab^ = t^\S'a + 5'it^\. (A.5) 

By a unitary change of basis we can bring into block-diagonal form and read off the 
irreducible representations. There are always two irreducible representations; the symmetric 
(antisymmetric) combination of two fundamentals has dimension Nd = Nc{Nc ± l)/2.^ 
Suppose we are interested in only one of these representations. We can assign it raised 
mid-alphabet latin indices, and from the change-of-basis matrix we can read off the Clebsch- 
Gordan coefficients K^ab- We can choose explicit values for the sextet of SU{3): 






(A.6) 



and for the antitriplet: 





1 

71 



/ 1 o\ 
-10 

y Oy 



(A.7) 



Note that, in the special case of SU{3), the antisymmetric combination of two fundamentals 
is an antifundamental; thus we can replace the raised index ' with a lowered index „ and write 
Kabc = ^abc/V^- It wiU also be useful to define Clebsch-Gordan coefficients for the conjugate 



^Whenever the symbols ± and =F are used in this Appendix, the upper (lower) sign refers to the symmetric 
(antisymmetric) combination of two fundamentals. 
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representation: Ki""^ = (K'^ba)* = K'^ha since we can always choose real coefficients. The 
coefficients K'^ab ^-I'e normalized such that 

K'K^ = Nd/Nc = {Nc ± l)/2, (A.8) 

and they satisfy the orthogonality relation 

TtK'Kj = 6i, (A.9) 

and the completeness relation 

K\bKi''' = liSi5t±5^jt). (A.IO) 
A. 3 Generators in the diquark representation 

We can form the generators in the diquark representation simply by extracting the ap- 
propriate block from the product representation: 

T'^'j = K\bK/'T^''\d. (A.ll) 

Applying Eq. ( [A.5 ) and simplifying, we find 

r^*j = 2 Tr KH'^Kj. (A. 12) 

Using this definition, it is straightforward to verify that the satisfy the appropriate com- 
mutation relation: 

[T^,T^]=if^^^Tc. (A.13) 
The quadratic Casimir operator in the diquark representation is given by 

T'^Ta = Cd = {Nc^I){Nc±2)/Nc. (A.14) 
The generators also satisfy the orthogonality relation 

TrT'^T^ = TdS^^ = ^{Nc ± 2)6^^, (A.15) 

and together with their products T^T^ they satisfy the completeness relation 



{T^T^f,{TATB)\ T ^^^T^',Ta\ - ^^^^^^ ^ ^i^f = S16^. (A.16) 

Due to its greater complexity, Eq. ( [A. 161) is not as useful as its counterparts Eqs. (|A.4D and 
( |A.10 ). Nevertheless, we include it here for the sake of completeness. 



- 16 - 



A. 4 Color factors 

Provided we sum or average over the color states of external particles, color matrices in 
the calculation of a cross section always appear with all indices contracted. Therefore it is 
straightforward to evaluate any color factor using Eqs. (^J), dOop , and ([A12|) . The color 
factors appearing in our calculation are: 

TrK^K. = Nn = ^^^i^, (A.17) 



{Nl-l){Nc± 2) 
2 

1 N?,-l 



T'^'.TAh Tr K'^K, = CdNd = ^ (A.19) 



Ti-KH'^K.itAf = ±-CfNc = ±^^, (A.20) 
T- Tri^^^.i^. ^ '-C.N, ^ (^'c-^fc±2) ^ ^^^^^^ 

B. Details of the NLO Calculation 

Here we give the details of the calculation of the NLO corrections to single diquark production. 
All divergences are regulated using dimensional regularization in 4 — 2e dimensions and the 
MS scheme is used to cancel the ultraviolet and collinear divergences. 

B.l Virtual corrections 

The Feynman diagrams contributing to the virtual NLO corrections are shown in Fig. [l|. The 
quark self-energy diagrams do not contribute since the quarks have no mass. Also, in the 
Feynman gauge, the contribution from the diquark self-energy diagram is found to be zero. 



Hence, as mentioned in Section 4.1, the self-energy diagrams and counterterms modify the 
Born cross section via the replacement |^2| 

A^ZA(Z|)-i(Zf)-V2A, (B.l) 

where Zl and are the wavefunction renormalization constants of the quark and diquark, 
respectively, Zx is the vertex renormalization constant, and the coupling constant on the 
right-hand side is renormalized. The expressions for the renormalization constants are given in 
Eqs. ( [4.1[) and (|4.2| ). All the following coupling constants are understood to be renormalized. 

The remaining virtual correction comes from the triangle diagrams. Performing the loop 
integrals, the invariant amplitude for the vertex corrections is 

a. 



A^vertex = -2V2i^iu' {p,)C^ {XrPr + XlPl) J^C, 

47r 



+{t\K^ +KHa)T-', -2-- 

1 Y4^y 



u(p2), (B.2) 
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where the color indices of the fundamental representation are suppressed and fiR is the renor- 
malization scale. Using the identities from Appendix the contribution to the cross section 
from vertex corrections is found to be 



a. 



•7 vertex \S) 



27r 



2Cf 



1 



1 + 



vr 



<7Born(s) 



is . 



(B.4) 



The sum of the Born and virtual diagrams is then 



O'Bomis) + O-vertex(s) 



'-6{1-t] 



z,{zl)-\z^)-y' 



+ K' 



(B.5) 



+2Cf { + ^ + llnAr + 1 



2e 



vr 



where ctq = 2-kNj:)}? /Nq. The renormalization constants have canceled the UV divergences, 
hence the remaining e poles are from coUinear and soft divergences. 

B.2 Real gluon emission 

The process 



q{pi) + q{p2)^g{k) + D{l), 



(B.6) 



as shown in Fig. also contributes to NLO corrections. The invariant amplitude for this 
process is 



Mqq = - 2V2i^l'I^gse'^{k)(|,iU^{pl)C\\RPR + \lPl) 



(B.7) 



{21 + ky 



m 



D 



2pi • k 



KHa 



2p2 • k 



r 



where the quark color indices have been suppressed and the diquark color indices have been 
left explicit. Squaring the amplitude and summing over color results in 



\M 



2t 



+ 1 



Cf- 



sm 



(B.8) 



where is the angle between the gluon and one of the initial-state quarks. 

The collinear and soft divergences of real emission are regulated by using the 4 — 2e 
dimensional Lorentz-invariant phase space: 



dPS% 



{27Tf-^'2Eg (27r)3-2^2/o 
(16^)-!+^ (1 -r)i^2e 

r(i - e) F 



5''-''{pi+P2-l-k) 



d0sm^-^'e 



(B.9) 
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To illustrate the cancellation of soft divergences and factorization of collinear divergences, 

(B.IO) 



the phase-space slicing method |23, U] is used. The soft region is defined by 



and the collinear region is defined by 

1 - |cos6l| < (5, (B.ll) 
where Eg = — r)/2 is the gluon energy, and Xmin and 5 are small arbitrary parameters. 

B.2.1 Soft region 

Using the upperbound on final state gluon energy in Eq. (IB.IOD , the contribution to the NLO 
cross section from soft gluon emission is 



<7soft 



1 1 



dPS'^\Mqg\''QCYXmin- Eg) 



s 2tt 



2Cf h 2 In^ Xmin - — 

e"^ e 6 



(B.12) 



+ Cd[-- 21nXmin + 2 



The single e poles in the diquark color dependent corrections and the double e poles cancel 
between the vertex corrections of Eq. (|B.5| ) and soft gluon radiation. Only the collinear 
divergences are left. 

B.2.2 Collinear region 

The collinear contribution to the cross section is then given by 



CTcol 



1 1 

2§m^ 



dPS'2\Mqq\^@{Eg - ^X^in)G(5 - 1 + | COS 0| 



asCp 



vr 



1 1 ^ 
- +ln- 

e 2 



1 ^'min 



dxaBornixs) 



1 + -e{l-xf 

(1 - x)l+2^ ' 



(B.13) 



where 1 — x is the energy fraction of the gluon with respect to an initial state quark and the 
lower bound Eg > -^Xmin is imposed to avoid double counting with the soft gluon emission 
contribution. 

Performing the x-integration, the contribution to the real-gluon emission cross section 
from the collinear region is 



<7col = Z — <^Born[S) 



2 In Xr, 



2 In - In Xmin - 2 In^ Xr_ 



(B.14) 



+ / dxaBovn{xs) 
Jo 



1 S\ 1 + x^ 

- + In - ^ 

6 2;(l-x)+ 



+(l-x) + 2(l + x^)fMl^ 
V 1 — X 
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where the "plus distribution" is defined by 



/' dxg{x) [F{x)\ = [' ^dx g{x)F{x) - g{l - /3) [' ^ dy F{y). (B.15) 
Jo Jo Jo 



Adding the contributions calculated so far one obtains 



CBorn + ^vertex + "^soft + ^col 



S 



<^(1 - ^) + JC-^I - ^ - §1^;^ - 1 + V - 21n^lnx^i„ ]S{1 - r) 



TT 
"3 



5_. 
2 



+ Cd( l-21na;^in- — )(5(1-t) 



2Cf 



1 5\ 1 + , , , o,/ln(l-r) 



(B.16) 



The soft and UV divergences have been cancelled, hence the remaining divergences are 
collinear. These collinear divergences can be absorbed into the definition of the pdfs. Using 
the MS scheme, the universal counterterm for the collinear singularities is 



(B.17) 



where ^u^? is the factorization scale and Pqqir) = Ci?[(l + r^)/(l — r)]+ is the DGLAP splitting 
function. The sum of all the contributions to the NLO cross section so far calculated is 



CT 



fBorn + f vertex + Csoft + <7col + ^col 



TT 



+ Cd 1 - 21na:;inin - 



2Cf y— - 1 - 21n - Inxmin + - In 
27r2 



3 ' -"'2 2"- 



+ 



agCp 



TT 



In 



2 4y(l-r) 



+ (l-r) + 2(l + r^ 



1 -T 



(B.18) 



The remaining dependence on the arbitrary parameters Xmin and 5 will be removed by calcu- 
lating the contribution from hard scattering. 



B.2.3 Hard scattering 



The remaining portion of the two-particle phase space, given by the constraints Eg > ^Xmin 
and 1 — I cos^l > S, corresponds to hard gluons in the final state. The contribution to the 
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cross section from hard scattering is 
1 1 



<7Hard 



2^ 



j dPS'^\ Mqq\''@{l - 5 - I COSe\)Q{Eg - ^X^in) 

l + x^ -e{l- x)2 



2CF\n--CD 



1 ^min 



dx <7Born(a^s) 



(1-x) 



l+2e 



(B.19) 



Adding all the cross sections together, the total NLO cross section with qq initial state is 



11 



(B.20) 



+2Cf 
-Cd 



l-|ln4M(l-^) + l-- 



B.3 Gluon initiated process 

Next-to-leading order QCD corrections also arise from the gluon-initiated process 

9iPi)+qiP2) ^ q{k)+D{l), 
shown in Fig. 0. The invariant amplitude for this process is 



(B.21) 



Mgg = -2V2ig,fiy^{pi)(f>iU^ip2)C^{XRPR + XlPl) 



i (2/-Pi) 

s 2k -pi' " 21 -pi 



tAK Y ^ 



-7^ - K^Ta 



u{k), (B.22) 



where, again, the diquark color indices are left explicit while the quark color indices are 
suppressed. Performing the spin sums the invariant amplitude squared is found to be 



+ 2Cd 



1 -e 

1 — cos 6 \ l — T 
4t 



1 



2t ) -(3 + cos6')(1-t) 
+ 



(B.23) 



(1 + r)(l + /3cos6l) (1 + t)2(1 + /3cosi 

where /3 = (1 — r)/(l + r), and 9 is the angle between the initial state gluon and the final 
state antiquark. The cross section for the gluon-initiated process is then 

1 



1 



a 



4(1 - e)Nc{N^ - 1) 2s 



\Mgg\'dPSI. 



(B.24) 



There are collinear divergences but no soft divergences in Eq. ( |B.23|) . Following the 
same method as for gluon emission, we isolate the collinear divergences using the cutoff in 
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Eq. ( |B.11| ) . The contribution to the gluon-initiated production cross section from the cohinear 
piece is given by 

5coi = 4(i_,)^^(^._i) j| / dPSl\M,,\MS - 1 + I cos^l) (B.25) 



s 47r 



fl-T)2+r2 



-+ln^ + 21n(l-r)-l 



Again, the colhnear singularity can be absorbed into the definition of the pdfs. Using the 
MS scheme, the universal collinear counterterm is 



~CT _ cro as / ^T^fJ-n 



2 \ e 



s 27r\ fil J r(l - e) e 



1 1 



(B.26) 



where Pgq{T) = ^[(1 — t)^ + r^] is the DGLAP splitting function. Adding the collinear piece 
to the counterterm, we have 



0-col + 0-col 



s 47r\ 



+ 



In 



6s 
2^4 



+ 21n(l -r) - 1 



(B.27) 



The hard scattering cross section does not contain any singularities and can be calculated 
in 4 dimensions. The contribution to the cross section from hard gluon-quark scattering is 
then 



^Hard 



s Att 



j dPS2\Mgq\''Q{l -5-\ cos9\) 



2^ ^ 2 



(B.28) 



+ 



Cd 
Cf 



(1 -r)(l + 2t) +2TlnT 



Adding all the contributions together, we obtain the giuon-initiated cross section 

{1-t){7t-3) 



cr. 



— — CT ^ 

gq — 0"col + O'col + '''Hard 



2N^s 



Cd 
Cf 



In 



m 



D_ 



+ 21n(l-T) 



+ 



(B.29) 



(l-T)(l + 2r) + 2rlnr 
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